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We investigated the spatial distribution of quantum fluctuations in a squeezed vacuum field, 
generated via polarization self-rotation (PSR) interaction of an ensemble of Rb atoms and a strong 
near-resonant linearly polarized laser field. We found that the noise suppression is greatly effected 
by the transverse profile of a spatial mask, placed in both the squeezed field and the local oscillator, 
as well as its position along the focused beam near the focal point. These observations indicate the 
spatial multi-mode structure of the squeezed vacuum field. We have developed a theoretical model 
that describes the generation of higher-order Laguerre-Gauss modes as a result of PSR light-atom 
interaction. The prediction of this model are in a good qualitative agreement with the experimental 
measurements. 

PACS numbers: 42.50.Lc, 42.50.Nn 


I. INTRODUCTION 

The ability to control the spatial mode composition 
of optical fields promises practical applications for many 
quantum technologies [1]. Eor example, multi-mode 
squeezing and entanglement can be used for quantum 
imaging [2, 3], quantum information multiplexing [4, 5], 
hyper-spectral encoding, etc. At the same time, the pres¬ 
ence of several spatial modes may be detrimental for 
quantum-enhanced measurements [6]. If the quadrature 
angle, corresponding to the maximum quantum noise 
suppression, is different for different spatial modes, the 
overall detected noise may be increased, even if each in¬ 
dividual spatial mode is squeezed. However, if the proper 
mode structure is identified, it may be possible to extract 
a single squeezed mode and eliminate others by means of 
an efficient spatial mode sorter [7] or to tailor the appara¬ 
tus [8] to further boost the signal-to-noise ratio (SNR). In 
this way, we would be able to produce single, maximally 
squeezed, higher-order modes upon demand. 

In this paper, we investigate the spatial properties of 
the resonant optical field after interaction with dense Rb 
vapor under the conditions for polarization self-rotation 
(PSR) [9-1 1]. It is known that in the case of a linearly 
polarized input field, the noise in the orthogonal polar¬ 
ization is modified, leading to the generation of squeezed 
vacuum [12-14]. Several experiments demonstrated up to 
3 dB of quantum noise suppression with this scheme [15- 
17]. It is still an open question why better quantum noise 
suppression cannot be experimentally achieved, in con¬ 
trast to the predictions of the theory [12, 18, 19]. Here 
we present a study of the spatial mode structure of the 
output optical fields in the PSR squeezing process, that 
gives strong evidence of its multi-mode nature. In par¬ 
ticular, we carried out a simple sorting of the spatial 
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modes, elucidating the mode structure of the pump and 
the squeezed field. Eurthermore, we developed an intu¬ 
itive theoretical description of the higher-order spatial 
modes that qualitatively agrees with the experiment. 


II. EXPERIMENTAL APPARATUS 

The experimental setup for squeezing generation, de¬ 
picted at the Eig. I, is similar to the one used in our 
previous experiments [20]. An extended cavity diode 
laser was tuned by approximately 100 MHz to the red 
of the 5^5 'i/ 2,E = 2 ^ 5^Pi/2,E' = 2 transition of 
^^Rb (A 795 nm). The laser output was spatially fil¬ 
tered by passing it through a single-mode-polarization- 
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FIG. I. (Golor online) Experimental setup. SMPM is a 
single-mode polarization-maintaining fiber, A/2 is half-wave 
plate, GP is Glan-laser polarizer, PBS is a polarizing beam 
splitter, PhR is a phase-retarding wave plate, and BPD is 
a balanced photodetector. In the experiment a mask was 
inserted either in the collimated beam (a) or in the focused 
beam (b). 
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maintaining (SMPM) fiber, linearly polarized using a 
Glan-laser polarizer (GP) and focused into a 7.5 cm long 
cylindrical Pyrex cell with isotropically enriched va¬ 
por and no buffer gas. The waist of the focused beam (di¬ 
ameter 100 /im at 1/e^ intensity level) was located 6.5 cm 
from the front of the cell. The cell was mounted inside 
three layer /i-metal magnetic shielding, and maintained 
at a constant temperature of 66° G, corresponding to a 
Rb density of 5.4 x 10^^ cm“^. 

The interaction of a strong laser field with the ^^Rb 
vapor modifies the vacuum fluctuations in the orthogo¬ 
nal polarization. In order to measure quadrature noise, 
we rotated the polarization of the two fields by 45° and 
mixed them on a polarization beam splitter (PBS) be¬ 
fore sending them to a balanced photodiode (BPD). We 
thus realized a homodyne detection in which the strong 
laser held served as the local oscillator (LO) [15, 17]. 
The relative phase (quadrature angle) of the two fields 
was adjusted by tilting the phase-retarding (PhR) plate 
(a birefringent quarter-wave plate with a crystal axis set 
parallel to the LO polarization). All the measurements 
reported below were performed at 1 MHz detection fre¬ 
quency with 100 kHz resolution and 30 Hz video band- 
widths. To calibrate the SQL (Standard Quantum Limit, 
i.e., shot-noise) level, we inserted an additional PBS into 
the beam after the cell, which reflects the squeezed held 
and passes only the LO held. We detected noise suppres¬ 
sion of 1.9 ±0.2 dB below the SQL level in the maximally 
squeezed quadrature and 11 dB of anti-squeezing in the 
orthogonal quadrature (see Fig. 2 at 100% transmission 
level), which is quite typical for these kinds of atomic 
squeezers. 


III. SPATIAL PROPERTIES OF THE 
COLLIMATED SQUEEZED BEAM 

All previous experimental and theoretical analysis of 
PSR squeezing assumed an identical single spatial mode 
for both the strong pump and vacuum fields, with either 
fundamental Gaussian [15, 17] or Laguerre-Gaussian [20] 
transverse profiles. However, as we show below, the spa¬ 
tial composition of the squeezed held mode is more com¬ 
plex. 

In order to investigate the spatial property of the 
squeezed held, we inserted different radially symmetric 
spatial masks into the collimated beam after the Rb cell 
(shown in Fig. 1(a)). We classify the beam masks as irises 
or disks, depending on whether they block the outer or 
inner part of the laser beam. We used an adjustable 
iris, and a set of fixed sizes disks. These masks were in¬ 
serted into the collimated part of the laser beam after 
interaction with the atoms, making sure that the masks 
are properly centered. Since the masks reduced the LO 
power, we carefully recalibrated the shot-noise level for 
every mask and adjusted the phase retarding plate to 
track the maximally squeezed and anti-squeezed quadra¬ 
ture noises. 



Transmission Transmission 

FIG. 2. (Color online) Measured minimum (diamonds) and 
maximum (circles) quadrature noise when the laser beam is 
partially blocked by (a) an iris mask, (b) a disk mask, (c,d) 
the combined masks, formed by both iris and the disk. For 
all measurements the horizontal axis indicates the fraction 
of LO intensity, transmitted through the mask. In (c) the 
central disk alone blocked 8% of the LO power and in (d) the 
power loss was 25%. Dashed and dotted lines indicate the 
prediction of the single-mode model, described by Eq. (1). 
The zero of the vertical axis corresponds to the shot-noise 
noise level . 

The modifications of the quantum noise by different 
types of masks are presented in Fig. 2. Based on a single¬ 
mode description of the optical held, one would expect 
the changes in the detected quantum noise in all cases to 
depend only on the total optical transmission T, and to 
be accurately described by the beam-splitter expression: 

= 10 ■ logioT ■ + (1 - T)], (1) 

where SqV[n,out are the quadrature noise measured in 
dB before and after the mask. Clearly, the experimen¬ 
tally measured noise values deviate significantly from the 
naive expected dependencies; shown in Fig. 2 as dashed 
and dotted lines, indicating non-trivial spatial correla¬ 
tions [21, 22]. For example, even small losses of about 
10% (for the iris mask) and 30% (for the disk mask) bring 
the squeezed quadrature noise significantly above shot 
noise. Moreover, for the disk mask, even at small trans¬ 
mission (T < 5%), the quantum noise in both quadra¬ 
tures is more than 5 dB above shot noise, see Fig. 2(b). 
In contrast, Eq. (1) predicts that in all cases we should 
expect the noise to approach the SQL monotonically from 
below, never exceeding the shot noise level. 

To gain additional insight about the spatial distribu¬ 
tion of the squeezed vacuum held, we looked at the noise 
of a ring-like slice of the laser beam. To do this we con¬ 
structed a mask consisting of a fixed size opaque disk 
and a variable size iris. Fig. 2(c,d) represents modifica¬ 
tions of transmitted quantum noise by such masks where 
the fixed disk sizes are characterized by 8% and 25% 
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blockage of LO power. Again, Fig. 2 (c,d) shows that we 
were not able to improve the measured noise suppres¬ 
sion below that of the unobstructed beam, even though 
the anti-squeezing noise followed the uniform loss model 
much better. Perhaps, such combined masks were able 
to block especially noisy spatial modes. 

The above observations suggest that the generated 
squeezed field consists of several spatial modes, some 
of which are “noisier” than others. These modes are 
expected to be radially symmetric due to the cylindri¬ 
cal symmetry of our setup. The higher-order Laguerre- 
Gaussian (LG) modes [23] are thus natural candidates to 
use for the decomposition of the LO and the squeezed 
beam. Such models are written as: 


on(^W[Z) 


Cl 


xLiq 


w{z) 

2r^ 

w{zy 


( 2 ) 


giZ0gi(2p+|Z| + l) 3.vct3.n{z/ zr) 


where I is the azimuthal index and p is the radial index 
for each mode, Ci^p = ^ 2 p!/ 7 r(|/| +p)! is a normalization 
constant, wq is the beam waist, w{z) = wq^/I + {zIzrY 

is the width of the optical field, L^p are the generalized 
Laguerre polynomials, zr = ttujq/X is the Rayleigh range, 
and k = 27r/A is the wave number. Because of the con¬ 
servation of the optical angular momentum, we expect 
that only the modes with I = 0 can be generated via the 
PSR process. 


IV. SEMI-CLASSICAL THEORY OF THE 
SPATIAL MULTI-MODE GENERATION 


Since the previous theory [ 12 , 18] was performed in 
the plane-wave approximation, we develop a model that 
accounts for the possibility of the higher-order spatial 
modes. Treating the light classically, we start with the 
inhomogeneous wave equation 


1 _ 1 


( 3 ) 


where P = ^(d) is the polarization induced in the 
medium, (d) = ^ k'mnpmn is the expectation of the 

dipole operator in an atomic level basis {|n)}, and 
are dipole moment matrix elements. We use the density 
matrix expression for the polarization, along with the 
slowly varying envelope and paraxial wave approxima¬ 
tions, to transform Eq. (3) into a propagation equation 
in terms of the envelope functions E and p\ 

. d i n . ~ ik N _ . . 

m,n 


Next we apply a simple model in which we treat the D 1 
line of ^^Rb as a double-A scheme [24]. We treat the in¬ 
put field as a superposition of two circularly polarized 


fields characterized by Rabi frequencies and U_, cor¬ 
responding to right- and left-circular polarizations. We 
solve for the density matrix elements, and since our pump 
field is in fact linearly polarized, we convert to the linear 
polarization basis and distinguish the propagation equa¬ 
tions for the two envelope functions in terms of the Rabi 
frequencies: 

= + ( 5 ) 


where L = (^ — ^V^), k is the wave number, k 
is the coupling constant, 70 is the decay rate, |U|^ = 
and A is the detuning. We further note 
that Qy is the Rabi frequency of the ^-polarized pump 
field and Clx is the Rabi frequency of the x-polarized vac¬ 
uum field. The homogeneous equation solved in cylindri¬ 
cal coordinates yields the LG family of solutions ui^p{r)^ 
given by Eq. ( 2 ). We obtain the mode structure of the 
output beam by avoiding a numerical calculation and 
proceeding with a weak scattering approximation in the 
following way. We assume on the right-hand side 
(r.h.s.) of Eq. (5) take the form of the input modes, i.e., 
1 ^ 1 ^ ^ ^yi and use the fact that 
<C Sy to simplify Eq. (5): 

el 

L^x,v = -uofi{r) [/t7o tTI uo,o{r)*uoft{r) 

IS I (6) 

4,y 

In this approximation we can regard the r.h.s. of Eq. ( 6 ) 
as the source of the Rabi frequency on the left-hand side 
(l.h.s.) of Eq. ( 6 ), and with this observation we define the 
appropriate sources px and py to simplify the notation, 
i.e., we transform Eq. (5) into LQx,y = Px,y{'k)- We can 
now use a Green function method [25] to write an integral 
expression for Clx^y'. 


nx,y = j r'dr'd(l>' K{r\r')VLl°^°{r') 

+ J dz' J r'dr'dfj)' G{f\f')px^y{r'), 


( 7 ) 


where represents the portion 

of the beam passing through the cell unaltered, 
is the propagator, and G{f\f') is the Green function for 
Eq. ( 6 ). An adjustable parameter, related to the ab¬ 
sorption and properties of the cell, characterizes the rel¬ 
ative strength of px,y to Thus, once A(r|r') and 

G{f\r') are known the problem is solved. The most effec¬ 
tive way to write the propagator (and consequently the 
Green function) is in terms of the LG modes. We define 
the propagator and Green function accordingly: 


I P 

G(r|r') = &{z — z')K{f\r'), 


( 8 ) 
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where 0 is the Heaviside step function. 

We further utilize the LG modes with an expansion of 
the sources px^y' 


Px,y{'^) — 

I p 

= J rdrdfj) ulp{f) Px,y{^- 


(9) 


We pause to point out that, since the propagation equa¬ 
tion inherits no (j) dependence from either the atoms or 
the input fields, the beam solutions will be limited to 
I = 0 modes. This is clearly seen through the integral 
for the ci^p coefficients by noting p^^y = px,y{T^z) and 
separating the (j) phase from the LG mode: 


/ p27T 

rdrpx,y{r,z)ui p{r,z) J dc/) (10) 


As one can see, the 0 integral vanishes for / 7 ^ 0. The 
beauty of this approach is that, as a consequence of the 
orthogonality of the LG modes, we are able to find an 
analytic solution which retains the mode structure of the 
field. Using (8—10) in (7) we arrive at the final solutions 
for the Rabi frequencies in the weak scattering limit: 


a 


x,y 


n 


homo 

x,y 


(f) 


/ dz'co,p{z'). ( 11 ) 
^ JceW 


One can perform these integrals with ease and our cal¬ 
culations show that the summation converges rapidly for 
p < 5. Now, we proceed by introducing a new experi¬ 
mental investigation which can expose the complicated 
LG phase structure predicted by the preceding theory, 
and hidden in the beam mode structure. 



Position of Iris (cm) Position of Iris (cm) 


FIG. 3. (Color online) Minimum quadrature noise vs. a 
variable iris position. The iris size was adjusted to maintain 
constant transmission of the LO beam for each trace as de¬ 
picted in the legend: 55% (a, circles), 70% (b, squares), 75% 
(c, crosses), and 95% (d, plus signs). 


V. SPATIAL PROPERTIES OF THE FOCUSED 
OUTPUT BEAM 


The phase factor (2p+ |^| + 1) arctan(^/ 2 (i^) of the gen¬ 
eral LG mode (Eq. ( 2 )) is referred to as the Gouy phase, 
and is responsible for a phase shift immediately around a 
beam focus. If our LO and squeezed field do in fact have 
a mixed LG structure, then we should observe asymme¬ 
tries about such a focal point; otherwise, the beam will 
have a simple Gaussian profile on both sides of the fo¬ 
cus. Referring back to Eq. (11), we see that the semi- 
classical theory predicts a mixture of 5 different Gouy 
phases. Therefore, we built a one-to-one telescope to 
create an extra focal point and displaced the iris mask 
down the telescope to investigate the LG structure. The 
modified setup is depicted in Eig. 1 (b). As we moved the 
iris along the beam, we change its size to maintain the 
same transmission of the LO beam, i.e., we kept the same 
ratio of the iris radius and the w{z). Therefore, our data 
is mostly effected by the different Gouy phase change in 
different spatial modes. We tracked the quantum noise 
in the maximally squeezed quadrature versus position of 
the iris. As one can see in Eig. 3, the quantum noise 
changes very drastically as the iris is moved right around 
the focal point. In the case that the iris size is large 
enough to allow high transmission, we see that the noise 
level does not change much as the iris is moved along 
the beam (see Eig. 3(d)). On the other hand, as trans¬ 
mission of the mask decreased, the noise in the squeezed 
quadrature went above shot noise (see Eig. 3 (a, b, and 
c)). Perhaps the most surprising is the noise dependence 
for the 75% transmission mask (see Eig. 3(c)), where one 
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FIG. 4. (Color online) Comparison of the experimental and 
theoretical (left and right column) dependencies of the LO 
beam transmission (top row) and the squeezed field minimal 
noise power (bottom row) on the iris position for several fixed 
iris sizes. The legend denotes the peak transmission for each 
iris. The Rayleigh range is Zr =2.5 cm. To calculate trans¬ 
mission and noise power, we used Eq. (11) and Eq. (1). 
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can see a very sharp drop of the noise into the global 
minimum around z=l cm, where it is squeezed below 
the shot-noise level. Similar but less pronounced behav¬ 
ior can be observed for the 70% transmissive mask (see 
Fig. 3(b)). 

For the same setup (see Fig. 1(b)), we mapped the 
transmission and minimum noise curve vs. the position 
of irises with fixed diameters, as shown in Fig 4. The 
transmission of LO exhibits some degree of asymmetry. 
The position of the maximum transmission shifted as we 
varied the iris size (see Fig. 4(a)). For the irises with 93% 
and 86% peak transmission, there are more than one lo¬ 
cal minima. We use Eq. (11) to predict the transmission 
level of the two LO field (see Fig. 4(b)) and it is in a qual¬ 
itatively good agreement with our experimental data (see 
Fig. 4(a)). We would like to again emphasize, that it is 
the interplay of Gouy phase shifts in the mode superpo¬ 
sition (Eq. (11)) that create this peculiar dependence on 
the iris position around the focal point. 

Based on the beam splitter model (Eq. (1)), we cal¬ 
culate the expected minimum noise as in Eig. 4(d). In 
the theoretical plot, the noise is always lower than shot 
noise since the model does not take into account the ex¬ 
cess noise. But, we can still compare it with Eig. 4(c) 
and find the overall behavior is quite similar. In both 
plots, the minimum in the noise power traces shift to the 
right, i.e., farther from the cell as the iris size shrinks. 
The qualitative agreement of these plots is further evi¬ 
dence of the LG structure of the squeezed vacuum beam. 
A more rigorous description necessarily requires a quan¬ 
tum treatment of both the light-matter interaction and 
the interaction with the spatial masks. Therefore, we 
next develop a simple second quantized theory that in¬ 
corporates the LG structure of the semi-classical beam 
and predicts the distribution of the quantum noise for 
the individual LG modes. 


VI. MULTI-MODE QUANTUM NOISE 
CALCULATIONS 


the LO, simulate the results of the homodyne detection. 

We proceed with the conventional single mode second 
quantization procedure, modifying it along the way to in¬ 
corporate the higher-order LG mode structure. It is com¬ 
monplace to represent the atomic polarization by a power 
expansion in the applied electric field strength E [26] : 

oo 

P = E + X^^^E^ + ...). ( 12 ) 

n=l 


Typically, in the plane wave approximation, one would 
define the order susceptibility as = P^'^^/eoE^ 
[23]. However, in the paraxial approximation, special 
care must be taken since the electric field is confined to 
the beam axis. We reduce the electric field and polar¬ 
ization to 2 ) dependence by integrating over a slice of the 
vapor cell perpendicular to the light beam. This is nec¬ 
essary, since the electric field dies off rapidly in r, and 
valid, since we are well within the validity of the paraxial 
approximation, i.e., each modes transverse profile sim¬ 
ply rescales with Thus, from this point forward we 
are careful to distinguish the strictly 2 )-dependent elec¬ 
tric field and polarization 


— / rdrd(j)P{r, 0, z) 

icell 

— / rdrd(j)E{r^(j)^ z). 
JceW 


(13) 


Now that we have carefully chosen this representation, 
we may define the order susceptibility as 




P(”)(z) 
E'^iz) ’ 


(14) 


where i^z) is the order polarization in response of 
the input field E{z). Then we further expand the optical 
polarization in terms of the LG modes: 


(r, 4>,z) = ^ ci^pUi^p, (15) 

hv 


We now extend our previous semi-classical treatment 
of the multi-mode light-atom interactions to develop a 
second quantized theory that predicts the underlying 
modal structure of experimentally realized squeezing. At 
the first glance, one can think that our homodyne de¬ 
tection scheme is immune to any spatial mode structure 
modifications since the strong field, used as the local os¬ 
cillator (LO), also passes through the vapor cell, so any 
spatial profile distortion are shared by both the squeezed 
optical field and the local oscillator. However, we still 
need to take into account the mode structure. Even if the 
overlap between squeezed modes and the local oscillator 
modes is perfect, distinct modes are going to contribute 
into the measured noise level differently. Therefore, to 
predict the resulting measured quantum noise, one must 
first predict the squeezing parameters for each of the vac¬ 
uum mode, and then, using the semi-classical solution for 


where ci^p = J rdrd(p Ui^ P^^\r^ 0, z). Motivated by this 
expansion, we introduce the accompanying LG suscepti¬ 
bility for each LG mode: 



where we define the LG polarization P^^^z) = 

/cell ci^pUi^p. In this representation, the LG suscep¬ 

tibility not only carries information about the LG spa¬ 
tial structure, it can be interpreted as the mechanism for 
“cross-talk” between the Gaussian pump and the vacuum 
LG modes during the non-linear interaction. Therefore, 
we define an interaction Hamiltonian density following 
the prescription of Ref. [27]: 

« ifixf^l{z){alpP^ - aj^P) 


(17) 
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where ai^p is the l^p spatial mode operator of the vacuum 
and b is the pump input mode. Thus, the full interaction 
Hamiltonian is 


Hoc 



(18) 


We proceed by assuming the ^-polarized pump mode is 
an undepleted coherent state and make the substitutions 
b —and b^ —This leads to the 
evolution operator 


U = exp 



(19) 


where iji^p = Next, making the substitution 

J0i,p ^20) 


^i^p = 2 dz r]i^pt = 
JceW 


we transform Eq. (19) into the form of the familiar 
“squeezing” operator with the addition of the LG in¬ 
dices [28]: 


5(6 ^ U = exp [^ 1(6, al - ag) 

hv 


( 21 ) 


An inspection of Eq.(20) reveals that the behavior of the 
LG modes through the focus in the vapor cell drastically 
effects the distribution of squeezing among the modes. 
Likewise, the Gouy phase of each mode and the cell posi¬ 
tion/dimensions effectively determines the squeezing an¬ 
gle of each mode. 

After defining the proper squeezing operator, we can 
find the variances of the LG quadrature operators: 

i,p = Q T 

1 . .t ■ 

The calculation is not particularly straightforward, but 
the result [29], 


((AXi,2 6p)^> 


1 

4 


cosh^ ri^p + sinh^ ri^p 


^ 2 sinh ri^p cosh ri^p cos Oi^p 


(23) 


is the familiar textbook result, but written for a particu¬ 
lar LG mode [28]. In other words, in this simplified sce¬ 
nario, each mode of squeezing has the statistics of a single 
squeezed mode. Eurthermore, the photocurrent (differ¬ 
ence) variance is a weighted sum of quadrature variances 
(some squeezed, some anti-squeezed) of all modes which 
overlap the LO [30]. As we have shown previously, only 
the modes with I = 0 are relevant to our experiment due 
to symmetry and consideration of angular momentum. 
Then one can show 

{{Aidf) oc y] cos^(Arg Op) 

P (24) 

+6^'"*’ sin^(Arg Op)], 


where we define the overlap integral Op as 

Q _ I 

V'/lfiLoPdv’ 


(25) 


Olo is the local oscillator, and (pp is the spatial function 
of the mode (in this case simply the UQ^p LG mode). 
To visualize these results we use the squeezed vacuum 
Wigner function; labeled here for each /,p mode. 


Wi,p{x, y) = - exp - ip + a;^) coshn,, 

^ ^ (26) 
+ ((a;^ - ) cos di ,, + 2xy sin 9i ,,) sinh r; 


Using quantum state tomography, we have recon¬ 
structed the Wigner function for our squeezed vacuum 
state (Eig. (5a)). This Wigner function corresponds to 
noise suppression of 1.9 ± 0.2 dB below the SQL level in 
the maximally squeezed quadrature (the best suppression 
we were able to achieve in our experimental geometry), 
a modest amount of squeezing compare to the predic¬ 
tions of the single-mode theory [12, 18, 19]. However, the 
preceding theory suggests that more significant squeez¬ 
ing could be hiding in the spatial mode structure of the 


W(x,y) 



X 


FIG. 5. Wigner function of (a) the experimentally realized 
squeezed state and (b) a hypothetical minimum uncertainty 
state with 1.9 dB of squeezing along the Xi quadrature, shown 
here in 3D and contour. The simple second quantized theory 
predicts the multi-mode structure (Fig. (6)) that results in 
the same level of measured squeezing via homodyne detec¬ 
tion. The axis labels x and y are proportional to the Xi 
and X 2 quadratures respectively. The Wigner function has 
been rescaled by the peak amplitude of the vacuum Wigner 
function. 
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beam. To facilitate the comparison with the experimen¬ 
tal results, we employ the second quantized theory to de¬ 
termine the value of the coupling constant k in Eq.5 so 
that the minimum quadrature noise of the theoretically 
calculated output measured via homodyne detection pro¬ 
vides 1.9 dB noise suppression in the Xi quadrature, i.e., 
0/2 = 0. The Wigner function of such a hypothetical 
state is plotted in Fig. (5b). The same simple second 
quantized theory also gives us the decomposition of this 
state into various LG modes: first, we follow the pre¬ 
ceding progression using the polarization derived in the 
semi-classical section (r.h.s. of Eq. (5)) and the solu¬ 
tion fly (Eq. (11)) for the y-polarized beam since it is 
used as the LO during detection. Then, we re-scale the 
coupling constant such that the weighted sum of quadra¬ 
ture variances (Eq. (24)) reproduces the same squeezing 
of our hypothetical state. To take care of the scaling, 
we first recognize that since 0 = 0, Eq. (23) reduces to 
((AXi^ 2 /,p)^) = (l/4)e^^’^^’. Furthermore, in the limit of 
a strong LO, the photo-current difference variance is 

((Aid)2)oc4|fLOp((AX(0))2) (27) 

where X{9) = (l/ 2 )(ae“®® + a’^e®^) is the field quadrature 
operator at the angle 0. Therefore, combining Eq. (24) 
and Eq. (27) for the special case ^ = 0, we impose that 
the scaling parameter 7 satisfy the equation 

I^LOpe-^'-exp = |Q^|2[g-27rp cos2(Arg Op) 

P (28) 

sin^(Arg Op)] 

where rexp ~ 0.565 is the experimental squeeze parame¬ 
ter and I^LOp = f Olo- We solve for 7 and 

determine the re-scaled LG squeeze parameters = 7 
and point out that the squeeze angles Op = Arg(^p) re¬ 
main the same. 


TABLE 1. Squeezing Parameters for Various Modes 


p 

r' 

' p 

Op 12 

\Op\ 

Arg(C>p) 

0 

1.297 

160° 

0.995 

71° 

1 

0.315 

113° 

0.091 

101° 

2 

0.149 

97° 

0.031 

123° 

3 

0.029 

25° 

0.006 

76° 

4 

0.011 

171° 

0.004 

CO 

00 

0 

5 

0.010 

18° 

0.002 

160° 


Figure ( 6 ) shows the Wigner functions for the squeezed 
LG modes, which comprise the hypothetical state, the 
p = 3 and higher modes are omitted since they show no 
appreciable squeezing. Table I summarizes the squeez¬ 
ing parameters and overlap integrals for different modes. 
It is important to point out that if only the fundamen¬ 
tal mode p = 0 could be isolated, its minimum quadra¬ 
ture noise would have been measured to be more that 
11 dB below the shot noise. The p = 1 by itself (if iso¬ 
lated) would have displayed —3 dB of squeezing, but at 
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FIG. 6. The Wigner functions for the p = 0 — 2 modes 
((a),(b),(c) respectively) which, when measured simultane¬ 
ously via homodyne detection, recreate the example Wigner 
function depicted in (Fig. (5). The p = 3 and greater modes 
are omitted since they appear essentially as vacuum modes. 
The axis labels x and y are proportional to the Xi and X 2 
quadratures respectively. The Wigner functions have been 
rescaled by the peak amplitude of the vacuum Wigner func¬ 
tion. 


the different quadrature angle, and the trend continues 
for higher modes, quickly approaching a coherent state. 
Thus, we find that there is actually much more squeez¬ 
ing available in the individual modes, but in combination 
the measurable noise suppression is much worse due to 
the fact that the squeezing angles are out of phase. This 
deterioration is independent on the atomic excess noise 
[18, 19], that is not taken into account in our simple the¬ 
oretical model and that can further reduce the measured 
squeezing value. We hypothesize that the different Gouy 
phase for each mode rotates the squeezing angle for that 
mode at a different rate. 





















VII. CONCLUSION 

We demonstrated that the PSR squeezer generates a 
multi-spatial-mode squeezed field. We developed a semi- 
classical model, which qualitatively describes our exper¬ 
imental observation. Additionally, we present a simple 
second quantization procedure that shows how quadra¬ 
ture noise can be varied between LG spatial modes. Our 
results have application to improved single-mode squeez¬ 


ing and the production of squeezed light in higher-order 
modes upon demand. 
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